ABSTRACT
In this work, we develop mathematical models of the immune response to respiratory viral infection, taking into account some particular properties of the SARS-CoV infections, cytokine storm and vaccination. Each model consists of a system of ordinary differential equations that describe the interactions of the virus, epithelial cells, immune cells, cytokines, and antibodies. Conventional analysis of the existence and stability of stationary points is completed by numerical simulations in order to study the dynamics of solutions. The behavior of the solutions is characterized by large peaks of virus concentration specific to acute respiratory viral infections. At the first stage, we study the innate immune response based on the protective properties of interferon secreted by virus-infected cells. Viral infection down-regulates interferon production. This competition can lead to the bistability of the system with different regimes of infection progression with high or low intensity. After that, we introduce the adaptive immune response with antigen-specific T- and B-lymphocytes. The resulting model shows how the incubation period and the maximal viral load depend on the initial viral load and the parameters of the immune response. In particular, an increase in the initial viral load leads to a shorter incubation period and higher maximal viral load. The model shows that a deficient production of antibodies leads to an increase in the incubation period and even higher maximum viral loads. In order to study the emergence and dynamics of cytokine storm, we consider proinflammatory cytokines produced by cells of the innate immune response. Depending on the parameters of the model, the system can remain in the normal inflammatory state specific for viral infections or, due to positive feedback between inflammation and immune cells, pass to cytokine storm characterized by the excessive production of proinflammatory cytokines. Finally, we study the production of antibodies due to vaccination. We determine the dose-response dependence and the optimal interval of vaccine dose. Assumptions of the model and obtained results correspond to the experimental and clinical data.
ABSTRACT
Replication of viruses in living tissues and cell cultures is a "number game" involving complex biological processes (cell infection, virus replication inside infected cell, cell death, viral degradation) as well as transport processes limiting virus spatial propagation. In epithelial tissues and immovable cell cultures, viral particles are basically transported via Brownian diffusion. Highly non-linear kinetics of viral replication combined with diffusion limitation lead to spatial propagation of infection as a moving front switching from zero to high local viral concentration, the behavior typical of spatially distributed excitable media. We propose a mathematical model of viral infection propagation in cell cultures and tissues under the diffusion limitation. The model is based on the reaction-diffusion equations describing the concentration of uninfected cells, exposed cells (infected but still not shedding the virus), virus-shedding cells, and free virus. We obtain the expressions for the viral replication number, which determines the condition for spatial infection progression, and for the final concentration of uninfected cells. We determine analytically the speed of spatial infection propagation and validate it numerically. We calibrate the model to recent experimental data on SARS-CoV-2 Delta and Omicron variant replication in human nasal epithelial cells. In the case of competition of two virus variants in the same cell culture, the variant with larger individual spreading speed wins the competition and eliminates another one. These results give new insights concerning the emergence of new variants and their spread in the population.